function empirical_cov_list = ecf_over_nu(hat_theta, rep_base_se2, hat_vartheta, rep_val_se2, meta_nu_list)
%INPUT: baseline and validation study data, list of meta analyst's nu
%OUTPUT: empirical coverage frequencies evaluated at each nu in meta_nu_list
    
    % Initialization
    K = size(meta_nu_list, 2);
    empirical_cov_list = zeros(1, K);
    
    for k = 1:length(meta_nu_list)
        % Set the value of meta analyst's nu
        meta_nu = meta_nu_list(k);
    
        % Compute posterior means and variances for tau
        [bar_tau, bar_V_tau] = compute_posterior_moments(hat_theta, rep_base_se2, meta_nu);
        
        % Check whether the validation study estimates are located in the predicted interval
        empirical_cov = check_CI_over_nu(hat_vartheta, rep_val_se2, bar_tau, bar_V_tau, meta_nu);
    
        empirical_cov_list(1, k) = empirical_cov;
    end
end